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Abstract 

We show that if a Levy process creeps then, as a function of u, the renewal 
function V(t,u) of the bivariate ascending ladder process (L~ l ,H) is absolutely 
continuous on [0, oo) and left differentiable on (0, oo), and the left derivative at u is 
proportional to the (improper) distribution function of the time at which the process 
creeps over level u, where the constant of proportionality is d^ 1 , the reciprocal of 
the (positive) drift of H. This yields the (missing) term due to creeping in the 
recent quintuple law of Doney and Kyprianou (2006). As an application, we derive 
a Laplace transform identity which generalises the second factorization identity We 
also relate Doney and Kyprianou's extension of Vigon's equation amicale inversee 
to creeping. Some results concerning the ladder process of X, including the second 
factorization identity, continue to hold for a general bivariate subordinator, and are 
given in this generality. 

Keywords: Levy process, quintuple law, creeping by time t, second factorization identity, 
bivariate subordinator 
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1 Introduction 

Let X = {X t : t > 0}, X = 0, be a real-valued Levy process with characteristic 
triplet (7,<T 2 ,ilx), thus the characteristic function of X is given by the Levy-Khintchine 
representation, Ee %eXt = e mx ^ s \ where 

tf x (0) = ifl 7 - a 2 9 2 /2 + [ (e Wx - 1 - Wxl { \ x \ <1} )IL x (dx), for 9 e R, t> 0. (1.1) 

Jr 

X is said to creep across a level u > if P(t u < oo, X Tu = u) > where 

t u = inf{t > : X t > u}. 
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Our initial interest is in the time at which X creeps. Thus we introduce the (improper) 
distribution function 

p(t, u) = P(t u <t,X Tu =u), u>0,t>0. (1.2) 

We prove certain regularity properties of p(t, u) which allow us to relate it to the renewal 
function of the bivariate ascending ladder process of X. This yields the missing term, due 
to creeping, in Doney and Kyprianou's (2006) quintuple law. Using the quintuple law, 
we derive a Laplace transform identity which generalises the second factorization identity 
due to Percheskii and Rogozin (1969). We also relate creeping to Doney and Kyprianou's 
extension of the equation amicale inversee of Vigon (2002). Some of these results ex- 
tend from the bivariate ladder process to general bivariate subordinators, and we develop 
several of the results in this setting. In particular, it appears to have gone previously 
unnoticed, that the second factorization identity is a special case of a general transform 
result for bivariate subordinators. The results in the fluctuation setting are stated in 
Section [3j with their proofs given in Sections [51 The general bivariate subordinator case 
is developed in Section HI 

By a compound Poisson process we will mean a Levy process with finite Levy measure, 
no Brownian component and zero drift. The indicator of an event A will be denoted by 
1a, or sometimes by 1{A), and we adopt the convention that the inf of the empty set is 
+oo. 



2 Fluctuation Setup 

We need some notation, which is very standard in the area. Let (L s ) s >o be the local time 
at the maximum, and (L7 1 , H s ) s > the weakfy0 ascending bivariate ladder process of X. 
Bertoin (1996), Chapter VI, and Kyprianou (2006), Chapter 6, give detailed discussions 
of these processes and their properties; see also Doney (2005). When X t — > — 00 a.s., 
(L _1 ,if) is defective and may be obtained from a nondefective bivariate subordinator 
(C^ l ,~H) by exponential killing at some rate q > 0, say. Thus 

{l:\h s ) = { { ^ 1 ' Us) iis<e ^ (2.1) 

I (00, 00) if s > e(q), 

where e(q) is independent of {C^.l-i) and has exponential distribution with parameter q. 
In the case that (£ _1 , H) is nondefective there is no need to introduce exponential killing 
and we set (L~\ H) = (C~ l , H). 

We denote the bivariate Levy measure of (£ _1 ,%) by Tl L -i H (-, ■), and its marginals 
by n L -i(-) and ILy(-). The Laplace exponent n(a,b) of (L~ l ,H) is given by 

e -K(a,b) = E ^-aL~'-bH 1 . j < ^ = e ~ q ^-aC^ 1 -b Hl 



lr rhe distinction between weak and strict only makes a difference when X is compound Poisson. 
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for values of a, b G K for which the expectation is finite. It may be written 

K (a,b) = q + d L -ia + d H b+ [ [ (l - e -»*- fa ) H L -i >H {dt, dx), (2.3) 

Jt>0 Jx>0 

where d^-i > and dn > are drift constants. The bivariate renewal function of [L _1 , H) 

poo poo 

V(t, x)= P(L; X <t,H s < x)ds = / e- qs P(C- s l <t,H s < x)ds. (2.4) 
Jo Jo 

It has Laplace transform 



e 



-at—bx 



t>0 Jx>0 



for all a, b such that n(a, b) > 0. The positivity condition on k clearly holds when a, b > 
and either aV6>0org>0. 

Let X t = —X t , t > denote the dual process, and {L" 1 , H) the corresponding strictly 
ascending bivariate ladder processes of X. This is the same as the weakly ascending 
process if X is not compound Poisson. The definition of (L _1 ,if) when X is compound 
Poisson is as the limit of the ascending bivariate ladder process of X t — et as e 4- 0. All 
quantities relating to X will be denoted in the obvious way; for example IIj-i _§(•,•)> 
«(•,•) and V(-, •). We choose the normalisation of the local times L and L so that the 
Weiner-Hopf factorisation takes the form 

n(a, 0)/c(a, 0) = a, a > 0. (2.6) 

This would not be possible in the compound Poisson case if (L -1 ,if) were the weak 
bivariate ladder process; see Section 6.4 of Kyprianou (2006). 

3 Creeping Time 

It is well known that X creeps across some u > iff X creeps across all u > 0, in 
which case we say that X creeps. A necessary and sufficient condition for creeping is that 
dn > 0; see Theorem VI. 19 of Bertoin (1996). Our first result describes the (improper) 
distribution function of the time at which X creeps across level u. 

Theorem 3.1 (Creeping Time) 

(i) The following are equivalent: 

pit, u) > for some t > 0, u > 0; (3.1) 



2 Throughout, we will write II^-i •) and V(-,-) with the time variable in first position, followed 
by the space variable. This is at variance with some established literature, but seems desirable for 
consistency. 
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p(t,u) > for all u > and all t sufficiently large (depending on u); (3.2) 

p(t,u) > for all t > and all u sufficiently small (depending on t); (3.3) 

> 0. (3.4) 

(ii) If dn > 0, then for every t > 0, V(t, 0) = 0, V(t, ■) is absolutely continuous on [0, oo) 
with a left continuous left derivative on (0, oo), and for each u G (0, oo) satisfies 

p(t,u) = d H -^V(t,u) (3.5) 

where d-jd-u denotes the left partial derivative in u. 

(Hi) If df/ > and X is not compound Poisson with positive drift, then V(t, ■) is differ- 
entiable and p(t,-) is continuous on (0, oo) for each t > 0, and p(-,u) is continuous on 
[0, oo) for each u > 0. 

Remark 3.1 (i) Unlike the creeping case, it is possible that X creeps over some u but 
not all u by a fixed time t > 0. This is illustrated in Examples 15. II and 15.21 Nevertheless, 
from Theorem 13.11 X creeps over some u > by some time t > iff X creeps, and 
we obtain the generalisation (13. 5p of Kesten and Neveu's formula for the probability of 
eventually creeping over u; see pp. 119-121 of Kesten (1969). 

(ii) It follows immediately from Theorem 13.11 that 

d H V{dt, du) = P{t u G dt, X Tu = u)du. (3.6) 

This formula has already been noted by Savov and Winkel (2010), p. 8, and attributed 
to Andreas Kyprianou. Conversely, from Savov and Winkel's observation (13.61) . it follows 
that dnV(t, du) has a density given by p(t, u) for a.e. u. This however gives no information 
about p(t,u) for a given level u. One of the main points of Theorem 13.11 is that p(t,u) is 
the left derivative of dnV(t,u) for every u > 0, t > 0. This is particularly relevant in the 
quintuple law below. 

(iii) In the case that X is a subordinator which creeps, the Laplace transform of the 
time at which it creeps over u is given by 

E(e- ar -,X Tu =u) = d x v a (u) (3.7) 

where v a is the bounded continuous density of the resolvent kernel 

POD 

V a (du) = / e~ at P(X t G du)dt; 



o 



see page 80 of Bertoin (1996). This in principle gives the distribution of the time at which 
X creeps over u. Indeed 



d I°° d I°° P 

v a (u) = — / e- at P(X t < u)dt = — / ae~ at dt / P(X S < u)ds. 
Jo du J J 



du 
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Hence if the derivative could be moved inside the integral, from ( 13. 7p we would obtain 

P(t u < t, X Tu = u )=d x ^-^ P(X S < u)ds. (3.8) 

Since t 

d H V(t,u) = d x [ P(X s <u)ds 
Jo 

when X is a subordinator, it follows from (I3.5p that 

P(r u < t, X Tu = u) = d x -^— [ P(X S < u)ds. 

Thus (13.81) is correct provided d/du is replaced by cL/cL-u. Conversely one can use ( 13 .5ft 
to give an alternative proof of ( 13. 7D . 

(iv) Theorem 13.11 is concerned with regularity of V(t, ■). Some information about 
regularity of V(-, u) may be gleaned from Theorem 5 of Alili and Chaumont (2001), from 
which it follows that V(-,u) is absolutely continuous for each u > provided is regular 
for both (— oo, 0) and (0, oo), for X. 

The quintuple law is a fluctuation identity, due to Doney and Kyprianou (2006), 
describing the joint distribution of five random variables associated with the first passage 
of X over a fixed level u > when X Tu > u. Using Theorem 13.11 we are able to account 
for the contribution due to creeping, that is the term when X Tu = u. To give the result 
set 

X t = sup X s and G t = sup{0 < s < t : X s = X s }. 

0<s<t 

The quintuple law concerns the following quantities: 

• First Passage Time Above Level u: t u = mf{t > : X t > u}; 

• Time of Last Maximum Before Passage: G Tu ^; 

• Overshoot Above Level u: X Tu — u; 

• Undershoot of Level u: u — X Tu _; 

• Undershoot of the Last Maximum Before Passage: u — X Tu _. 

Theorem 3.2 (Quintuple Law with Creeping) Fix u > 0; then for x > 0, v > 0, 

< y < u A v, s > and t>0 

P (X Tu - u E dx, u - X Tu _ G dv , u — X Tu _ E dy, r u - G Tu — E ds, G Tu - E dt) 

d- 

= l {x>0} \V(dt, u - dy)\V(ds,dv - y)T[ x {dx + v) + d#- — V(dt, u)8 (ds, dx, dv, dy), 

O-U 

(3-9) 

where 5q is a point mass at the origin, and with the convention that the term containing 
the differential d_V(dt,u)/d_u is absent when dn = (in which case d_V(t,u)/d_u need 
not be defined). 
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The contribution to (13. 9 j) for x > is Doney and Kyprianou's quintuple law. Theorem 
13.21 then follows easily for a.e. u from Savov and Winkel's observation fl 3 . 6 j) . but this is 
clearly unsatisfactory, since it says nothing about any given u. To get the result for 
every u, (13.51) is needed. As a simple consequence of Theorem I3.2[ we record the joint 
distribution of the first passage time and overshoot of a level u > 0; 

Corollary 3.1 Fix u > 0. Then for x, r > 

P(X Tu - u e dx, t u e dr) = I(x > 0) / / |V(ds, « - d2/)|II L -i )H (dr - s,y + dx) 

J0<s<r J0<y<u 
8- 

+ d H - — V(dr,u)5 {0} (dx). 
In particular the distribution of the first passage time is 

P(r u edr)= [ [ \V(ds,u-dy)\U L -i >H (dr-s,(y,oo)) + d H ^—V(dr,u). 

J0<s<r J0<y<u c, - u 

Using the quintuple law, Doney and Kyprianou (2006) (Corollary 6) obtain the fol- 
lowing useful extension of the equation amicale inversee of Vigon (2002); for s > and 
x > 0, 

IT L -i^(ds,dx) = / V(ds,dv)Tlx(dx + v). (3.10) 

Jv>0 

They state this result for s > 0,x > 0, but their proof works equally well when s = 0. 
Observe that ( 13.101) gives no information about U L -i H on {(s, 0) : s > 0}. When X is not 
compound Poisson, U L ~i H (ds, {0}) is seen to relate to creeping as the following result 
shows: 

Theorem 3.3 Assume X is not compound Poisson. Then X creeps iff n^-i ^(ds, {0}) 
is not the zero measure. 

Despite the connection with creeping, it is easily seen that the jumps of (L~ l ,H) for 
which AH = do not occur when X creeps over a fixed level; see (I5.3P in Section |5j As 
a consequence of Theorem 13.31 we are able to characterise when ( 13. 10[) holds for all s > 
and x > 0: 

Theorem 3.4 (13. 10ft holds for all s > 0, x > iff X does not creep. 

When X is compound Poisson it does not creep, and we can deduce from Theorem 
S3] that 

n L -i jH (d S ,{0})= / V(ds,dv)Tl x ({v}), s>0. (3.11) 

Jv>0 

If IIx is diffuse then U L -i H (ds, {0}) reduces to the zero measure, but in general it may 
have positive mass. Thus Theorem 13.31 cannot be extended to the compound Poisson case. 
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The next result is an application of Theorem 13.21 to computing a quadruple Laplace 
transform. The finiteness conditions on k, below, clearly hold when p, p, A, u, 9 > 0, and 
in that case k,(u, p) > except when v = p = and q = in (12.31) . More generally 
the conditions allow for distributions with exponential moments, which can arise quite 
frequently in applications. 

Theorem 3.5 (A Laplace Transform Identity) 

Fix p, p, A, z/, 9 so that k(9, p + A), k(9, p) are finite and k(z/, p) > 0. 
(i) If A ^ p — p, then 



k(9, p + A) — k(9, p) 

(p + \-p)k{v,h) 

(3.12) 



f e'^E { e -p(Xr u -u)-\(u-x Tu _)-vG Tu —e(T U ~G Tu -). Tu < \ du 

Ju>0 ^ ' 

(ii) If A = p — p then 

[ E ( e -^-^-^,--«(r.-^-) ; Tu < oo) dn = , l -?±^M (3.13) 

A>0 V J K{u,fi) d+p 

provided the right derivative exists. 

The right derivative in (I3.13P exists and equals the derivative if k(9, p — e) < oo 
for some e > 0. When A = 0, 9 = v > 0, p > 0, p > 0, (13. 12[) reduces to the second 
factorization identity, (3.2) of Percheskii and Rogozin (1969). [Alili and Kyprianou (2005) 
give a short and elegant proof of the second factorization identity using the strong Markov 
property] Theorem 13.51 can be used in the computation of certain exponential Gerber- 
Shiu functionals from insurance risk, see Griffin and Mailer (2010a). Another application, 
to stability of the exit time, can be found in Griffin and Mailer (2010b). 

It is natural to ask if there is a quintuple Laplace transform identity, analogous to the 
quintuple law. It is straightforward to follow calculations similar to those in the proof of 
Theorem 13.51 and derive a corresponding expression to ( 13. 12ft . but because the component 
X Tu _ cannot be expressed in terms of the ladder process, the resulting expression cannot 
be expressed simply in terms of the kappa functions. 



4 Bivariate Subordinators 

(L _1 , H) and (£ _1 , H) are, possibly killed, bivariate subordinators, and some of our results 
require only this property. In this section we prove several theorems in this generality. 
These will then be applied in Section [5] to the fluctuation variables. 
In the fluctuation setting, let 

T u = = inf{s > : U s > u}, u>0. (4.1) 

Using H s = X L -i on {L^ 1 < oo}, s > 0, and recalling the exponential killing described 
in (12. ip . we have 

X Tu = H Tu} X Tu _ = H Tu -, r u = and G r „_ = on {T u < e(q)}. (4.2) 
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Thus, via ffL2|) . 

p{t,u) = P{C£ < t,Ur* =u,T u < e(q)). (4.3) 

This suggests the following setup. Let (Z, Y) be any two dimensional subordinator ob- 
tained from a true subordinator (Z, y) by exponential killing at rate q > 0, say. Corre- 
sponding to (12.31) . (Z,Y) has Laplace exponent n ZjY (a,b) = q — \ogEe~ aZl ~ byi where 



«z,y(a, b) = q + d z a + d Y b + [ [ (l - e~ at - bx ) U z , Y (dt, dx), 

Jt>0 Jx>0 



(4.4) 



for values of a, b G K for which the expression is finite. Analogous to the fluctuation 
variables, we define 

= inf{s > : y s > it}, u > 0, 

and 

p z , Y (t,u) = P(Z T y < t,y T y = u,T* < e(q)), (4.5) 
where e(q) is an independent exponential random variable with parameter q. Also set 

poo 

V z>Y {t, u) = / e-i s P(Z s <t,y s < u)ds. 
Jo 

So Vz )Y (-, •) has Laplace transform 

/ / e- at - bx V z , Y (dt,dx) = 1 -— if K ZtY (a,b) >0. (4.6) 

Jt>oJx>o K ZtY (a,b) 

We begin by investigating aspects of the regularity of p ZiY defined in 04. 5p . 

Lemma 4.1 The function p ZtY (-, •) has the following properties: 

(a) p z , Y (-,u) is right continuous and non-decreasing on [0, oo) for every u > 0; 

(b) p z>Y (t, •) is left continuous on (0, oo) for every t > 0; 

(c) If p z , Y {-,u) is continuous on (0, oo) for every u > 0, then p z , Y (t, ■) is continuous on 
(0, oo) for every t > 0. 



Proof of Lemma I4.lt First observe that the results trivially hold if d Y = 0, since then 
y does not creep and so p Z yit, u) < P(y T y = u) = for allt > 0, u > 0. Thus for the 
remainder of the proof we assume d Y > 0. 

Part (a) follows immediately from the definition of p ZjY . To prove Parts (b) and (c) we 
will use the following two equations which are simple consequences of the strong Markov 
property (cf. Andrew (2006)): for any x > 0, y > 0, r > 0, s > 0, we have 

Pz,v{r + s,x + y) > p z , Y (r, x)p z , Y (s, y) (4.7) 

and 

Pz,y{s, x + y) < p Z;Y {r, x)p ZjY {s, y) + l- Pz,v{r, x). (4.8) 
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By Theorem III. 5 of Bertoin (1996), which applies to nondefective subordinators with 
dy > 0, we have ]xai e yoP(y T y — s) = 1. Since dy > 0, y is strictly increasing, and so 
1 a.s. as e 1 0. Thus for every 8 > 

Iimp z , y («,e) = limP(Z ry < 5,y T y = e,T? < e(q)) = 1. (4.9) 
Now fix it > and £ > 0. Then for any < £ < it and 5 > we have by (14. 7p and 

Pz,y(*, w) - 1 + Pz,y(8, e) < Pz,y(£, e)Pz,y{t, u-e)< p z ,y{t + 5,u). 

Letting e J, then 5 J, 0, and using Part (a) and (14. 5J) , proves Part (b). Similarly if in 
addition, £ > and 6 <t, then 

p Zi y(5, e)p z ,y{t -6,u)< pz,y{t, u + e) < p z ,y{5, e)p z ,y{t, u) + 1 - P2,y(<5, e). 

Letting e 1 then 5 1 0, we conclude that 

Pz,y(t-,u) < liminf p z<Y (t,u + e) < limsupp Zj y(£, it + e) <p Z y{t,u). 

Thus if pz,y (•, u) is continuous on (0, oo) for every it > 0, then p z ,y(t, •) is right continuous 
on (0, oo) for every £ > 0. Combining this with Part (b) proves p z ,y{t, •) is continuous on 
(0, oo) for every £ > 0. □ 

Lemma 4.2 For any u > and £ > 0, 

p z ,y(£, v)du = d Y V zx (t, it). (4.10) 



Proof of Lemma 14. 2t If dy = then y does not creep, so both sides of (I4.10p are zero. 
Thus we may assume dy > 0. First observe that for any s, {v : y T y = v,T^ = s} is at 
most a singleton, and so 

POO 

/ 1{3^ = v, = s}dv = 0. (4.11) 

Next, if Tf = inf {s : Z s > £}, then 

{s < if} C {Z s < £} C {s < Tf}. (4.12) 
Thus, using KTT\f and fl4TT2|) . 

pit pU 

/ Pz,y(*, ^)dw = / P(2 T y < t,y T y = v,T^ < e(q))dv 
Jo Jo 

= P P(T* < if, y T y = v, T* < e(q))dv 
Jo 

= [ U p(y T y=v,Ty<T?Ae(q))dv. 
Jo 



Now since dy > 0, y is strictly increasing. Thus if y hits v then it does so at time T^. 
Hence 

f 1(3^0. = v, T-f < if A e(q)}dv = !% ATf Ae(?) - J] A?. 

° s<T*A3f Ae(<z) 

as each quantity represents the Lebesgue measure of the set of points in [0,u] hit by y 
by time if A e(g). Since y r = dyr + J2 s <r A34, this gives 

/ p ZtY (t,v)dv = dyE(T* ATf Ae(q)). (4.13) 
Jo 

But by f)4.12p (which also applies to y and T^) 

poo poo 

/ l{Z s <t,y s <u,s<e(q)}ds< / l{s < if, s < T%, s < e(q)}ds 
Jo Jo 

POO 

= l{s< Tf, s < 7?, s < e(q)}ds (4.14) 
Jo 

POO 

< / l{Z s <t,y s <u,s <e(q)}ds. 
Jo 

Thus by (Q3I) and (OIL 

/ ^)df = d y / P(Z S < t, 3^ s < w, s < e(q))ds = d Y V zy {t, u). 

Jo Jo 

□ 

Theorem 4.1 Parts (i) and (ii) of Theorem IJ.il /ioW precisely as stated with pzy in 
place of p, dy in place of du, and Vz.y in place ofV. 

Proof of Theorem I4.lt Since Vz, Y (t, u) > for allt > and u > by right continuity 
of (Z, Y), we have by fl4.10p and Lemma OJ (b) that 

dy > iff pz,y(t, u) > for some t > 0, u > 0. 

On the other hand 

dy > iff < -P^t " = u i Tu < e (<?)) = li m Pzy(t : u) for every u > 0. 

Combined with monotonicity of p{-,u) for it > 0, these give the equivalence of the subor- 
dinator versions of (13. ip . (13. 2p and (13. 4p . To complete the proof of Part (i) observe that 
the subordinator version of (13.41) implying the subordinator version of (I3.3P was proved in 
(14.91) . while the subordinator version of (13. 3p implying the subordinator version of (13. ip 
is trivial. 

If dy > then y is not compound Poisson, so Vz,y(i, 0) = 0. The remainder of part 
(ii) follows immediately from (I4.10p and Lemma [4.1( b). □ 
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Lemma 4.3 For every u > 0, 

P(AZ T y > o, Ay T y = 0) = 0. 



(4.15) 



Proof of Lemma 14. 3t If y is compound Poisson, then P{Ay T y = 0) = 0, so the result is 
trivial. Thus we may assume that y is not compound Poisson, in which case y is strictly 
increasing. Thus by the compensation formula, p7 of Bertoin (1996), for every e > 

P(Ay T y = 0, AZ T y >e) = Ej2 Hyt- = M > Ay t = 0, AZ t > e} 

t>o 

POO 

= n^,y((e,oo) x {0}) / P(y t = u)dt. 

Jo 

The last expression is by Proposition 1.15 of Bertoin (1996). □ 

The next result is a "quadruple law" for (Z,Y), in a similar spirit to the quintuple 
law. 

Theorem 4.2 (Quadruple Law) For u > 0, x > 0, < y < u, s > 0, t > 0, we have 
P {y T y — u E dx, u — y T y_ G dy, AZ T y E ds, Z T y_ E dt; < e(g)) 
= l {x>0} \V zx (dt,u - dy)\Il zx (ds,dx + y) + d Y - — Vz,y(dt, u)5 (ds, dx, dy), 

(4.16) 

with the convention that the term containing the differential d_Vz t y(dt, u)/d_u is absent 
when dy = (in which case dJVzy(t,u)jd_u need not be defined). 

Proof of Theorem 14. 2t Fix u > 0. By the compensation formula, we get for x > 0, 

< y < u, s > 0, t > 0, 

P{yT^ ~ u e dx, u — y T y^_ E dy, AZ T y E ds, Z T y_ E dt; < e(q)) 

= P{Ay T y Edx + y,u- y T y_ E dy, AZ T y E ds, Z T y_ E dt; Tl < e(q)) 
= E l{Ay r Edx + y,u- y r - E dy, AZ r E ds, Z T - E dt,r < e(q)} 



r>0 



poo 

= E l{u-y r - E dy,Z r _ E dt,r < e(q)}drIl Z:H (ds,dx + y) 
Jo 

= \V z , H (dt,u - dy)\IL ZjH (ds,dx + y). 

(4.17) 

Thus we have left to consider the case x = 0. 

First observe that from Proposition III. 2 of Bertoin (1996), it follows that 

P(y T y- < u = y-tf) = °> for a11 u > °- ( 4 - 18 ) 
11 



Now suppose dy > 0. By part (ii) of Theorem 14.11 

p z , Y (t, u) = d Y —V z , Y (t, u). (4.19) 

This, together with (14.15ft and ( 14. 18ft . shows 

P^t* = M ' n — ^Tu— e ^Z T y E ds, Z T y_ e dt; < e(g)) 

= P{y T y = u, Z T y e dt; T% < e(q))8 (ds, dy) ^ 2Q , 

= dy- — V zx (dt,u)5 Q (ds,dy) 

for < y < u, s > and t > 0. When dy = 0, ( 14.20ft continues to hold since the lefthand 
side of ( I4.20P is because y does not creep. Adding ( 14.20ft to ( 14.17ft then gives ( 14.16ft . 
□ 

The next result is a generalisation of Theorem 13.51 to the subordinator setup. The 
conditions on k z ,y are analogous to those on k in Theorem 13.51 

Theorem 4.3 (A Laplace Transform Identity for Subordinators) 

Fix p, p, A, v, 9 so that Kz >Y {6, t 1 + k z,y(@, p) ar ^ finite and K Z y{y, p) > 0. 

(i) If A 7^ p — p then 

e'^E f e ~ p{y T^ u) ^ u - y Ty-^ z T y-- SA \y. T y < e(g) ] du 

u>0 ^ ' 

_ kz,y(0,V + A) - k Z:Y (0,p) 
(p + A- p)k z ,y{v,v) 

(ii) If A = p — fi then 



(4.21) 



/ E ( e - pAy ^-^--^-- 9AZ ^ ; 2? < e(qj) du = ) ?±H^lfl (4.22) 

Ju>o V ' «z,y(f,A*) d+p 



provided the right derivative exists. 



Proof of Theorem 14.31 Taking the expectation over the set {y T y > u}, from ( 14.16ft we 
find 



it 



[ [ [ [ e-^-^lVz^u-dy^Uz^dx + y) u 23 ) 

J0<y<u Jx>0 Js>0 Jt>0 V ' 

[ [ If e- p{x - u+w) - x{u - w) - ut - es V ZtY {dt,dw)Tl Z y{ds,dx). 

J0<w<u Jx>u-w Js>0 Jt>0 
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Now take the Laplace transform of both sides of (14.231) . For A 7^ p — p, we obtain 

e _/m dii / / e -p(x-u+w)-\(u-w) 

u>0 Jo<w<uJx>u—w 

e -(n+\- P )u e -(p-\)w e -px du ^ 424 j 



ui>0 J x>0 J w<u<w+x 



w>0 Jx>0 



p — p — A 

Since we may clearly also include x = in the last integral, we then have 

Ju>0 ^ ' 

— T / / e-^^,y(dt,dw;) 

/ f e ~ es {e~^ +x)x -e- px )U z<Y (ds,dx) 

Js>0 Jx>0 



p- p 



X 

s>0 Jx>0 

k ZjY (0, p + A) - K^.y (#, p) - (// + A - p)d y 

(/i + A - p)K Z! y(v,p) 
Kz,y(d,p + A) - k z ,y(9,p) dy 



(// + A - p)k z ,y{v,h) K Z ,Y(v,fl) 

(4.25) 

by @3D and gSJ. 

If dy > then we need to add in the second term in (14.161) due to creeping. From 
(14.161) and part (ii) of Theorem 14.11 we have 

Ju>0 ^ ' 

= d Y [ e-» u [ e-^^-Vz^^du 

Ju>0 Jt>0 c) - u 

= d Y [ e-» u I e- ut V zx (dt,du) 

Ju>0 Jt>0 
dy 



kz,y{v, p)' 



Added to (|£25j) . this gives (1QT]) . 

Now consider the case where A = p — p. Let e > and set A' = p — p + e. Then 
kz,y(9,P + A') is finite since kz,y(@,p) is finite. Thus by (14.211) 

/ £ ( e -^-^--("-^ J-"^-^. T y < e(g) \ dn 

__ K ZtY (0,P + g) - kz,y(9,p) 
ek Z)Y (v, p) 
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Letting e \, and using monotone convergence completes the proof of (14.221) . 



□ 



Results similar to (I4.2ip can be found in Winkel (2005). Winkel's interest in bivariate 
subordinators is in modelling electronic foreign exchange markets and he does not make 
the connection with the ladder height process. As mentioned in the introduction, it 
appears to have gone previously unnoticed, that the second factorization identity is a 
special case of a general transform result for bivariate subordinators. 



5 Proofs for Section [3 

We now turn to the proofs of the fluctuation results from Section [3J Recall the definitions 
of p(t, u) and T u in (II. 2p and (14. ip respectively, and from (14. 3p that 

p(t,u) = P{C T l u < t,H Tu = U,T U < e{q)). 

In view of the correspondences (L _1 ,iJ) -H- (Z,Y), p(t,u) -H- pzy{t,u), and V(t, u) -H- 
Vz,y{t-, u), we can carry a number of results directly across from Section HJ In particular, 
the results of Lemma 14.11 and Lemma 14.21 hold with p in place of pzy and V in place of 

V Z ,Y- 

Proof of Theorem I3.lt Parts (i) and (ii) follow immediately from Theorem 14.11 For 
Part (iii), supposing X is not compound Poisson with positive drift, by Theorem 27.4 of 
Sato (1999), 

P(t u = t, X Tu =u)< P{X t = u) = 

for all u > 0. Consequently p(-,u) is continuous on [0, oo) for every u > 0, and hence by 
Lemma H~T| p(t, ■) is continuous on (0, oo) for every t > 0. Since p(0, •) = on (0, oo) this 
continues to hold for t = 0. By (I4.10p . this then implies differentiability of V(t, ■). □ 



Example 5.1 Let X t = t+Y t where Y is compound Poisson with Levy measure Hy (dx) = 
5{iy(dx) + 5{_i}(dx). Since X t = t + J2 s <t an< ^ ^2s<t ^ s an integer, we have that 
for any t E (0, 1), u > 

p(t, u) > iff u E U~ =0 (n, n + t]. 

Furthermore p(t,u) > e~ 2tu for u G (0,£]. Thus neither p(t,-) nor p(-,u) is continuous, 
and, using (I4.10p . V(t,-) is not differentiable. Finally, in contrast to monotonicity of 
p(-, u), p{t, ■) is not monotone. 

Example 5.2 Set X t = t — Y t , where Y is a subordinator and Hy (R) = oo. Then clearly 
p(t, u) = for u > t. Thus there is no hope of proving p(t, u) > for all t, u > even in 
the situation of (iii) of Theorem 13.11 

Before turning to the proof of Theorem 13.21 we need a preliminary result, generalising 
(I4.18p . which is surely well known, but for which we can not find an exact reference. 
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Lemma 5.1 Fix u > 0. 

(i) If X is not compound Poisson, then 

P(X t _ ^u,X t = u for some t > 0) = P(X t _ = u,X t ^u for some t > 0) = 0; (5.1) 
^m) For any £euy process X and u > 

P(X Tu _ < it, X Tti = u, r u < oo) = 0. (5.2) 

Proof of Lemma I5.lt (i) Use the compensation formula to write for u > 
P(X t _ ^ «, X = it for some t > 0) < E ^ H x t- + u, X t - + AX t = u} 

fOO 



POO P 

E dt 1{X_ ^ u, X t _ + £ = «}n x (eZ0 

/* /*OG 

/ / p(x_= M -0dtn x «). 



The last expression is when X is not compound Poisson, since the potential measure of 
X is diffuse by Proposition 1.15 of Bertoin (1996). Similarly for any e > 

P(X t _ = u, |X - w| > e for some t > 0) < E l{X t _ = u, \AX t \ > e} 

POO 

= U x ([-e,s] c ) / P(X t _=u) dt = 0. 



Letting £ — > completes the proof. 

(ii) Clearly we may assume P(X Tu = u) > else there is nothing to prove. In that 
case X is not compound Poisson. Since 

{X Tu _ < u, X Tu = u,r u < oo} C {X t _ 7^ u, X t = u for some t > 0}, 

f[5T2l follows from ([EI]). □ 

Proof of Theorem 13.21 Fix v > 0, < y < u A v, s > and t > 0. For x > this is 
Doney and Kyprianou's quintuple law. If x = 0, then from (14. 2p . (14.161) and (15. 2p 

P (X Tu - u G dx, u - X Tu _ G du, u - X Tu _ G dy, r u - G Tu - G ds, G Tu _ G dt) 

= P (H Tu - uedx,u- H Tu - G dy, A£ T J G ds, G dt; T u < e(qj) 5 (dv) 

= d H - — V(dt, u)5 (ds, dx, dv, dy). 

O-U 

□ 

Proof of Corollary I3.lt This follows easily from Theorem 13 . 2 1 and (I3.10p . Alternatively 
use Theorem 14.21 □ 
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Proof of Theorem 13. 3t Assume X is not compound Poisson. We can first easily dispense 
with the case that is irregular for (0, oo); because then, is irregular for [0, oo), so by 
construction, L _1 and H jump at the same times (see p. 24 of Doney (2005)). Hence 
U L -i H (ds, {0}) is the zero measure. On the other hand, by the strong Markov property 
at time r u , X does not creep over any u > 0. Thus the result holds in this case. 

We now assume that is regular for (0, oo), in which case the closure of the zero set of 
X — X is a perfect nowhere dense set with probability one; see for example the discussion 
on p. 104 of Bertoin (1996). Fix r 6 Q, r > and set Y t r = X r+t - X r . Let 



and 

By independence 



r r = inf{t > : Y t r > X r - X r } 
A r = {XjTr = X r — X r , X r — X r > 0}. 



P(Ar) = [ P{X Tu = u)P{X T — X r £ du). 

Ju>0 

Thus P(A r ) > iff X creeps. On the other hand, IT^-i t n(ds, {0}) is the zero measure 
iff P(AL~ l > 0, AH S = for some s > 0) = 0. The result then follows from the key 
observation that 

{AL; 1 > 0, AH S = for some s > 0} =' (J A r . 

To see this, assume s is such that AH S = and AL^ 1 > 0. Let r G (Ljl, L^ 1 ) fl Q. Then 
X r — X r > 0, X L -i = X r and X t < X r for Ljl < t < Lj 1 . Since a.s. the zero set of 
X — X contains no isolated points, it follows that off this P-null set YJ r = X r — X r . Thus 
A r occurs. Conversely fix r e Q and assume that Y^ r = X r — X r > 0. Choose s > 
so that r G (L~l, Lj 1 ). Then on A r , H s = X L -i = X r = H s _. Hence AH S = and 
AL; 1 > 0. □ 

As remarked earlier, the jumps of (L~ l , H) for which AH = 0, do not occur when X 
creeps over a fixed level. This follows from an application of Lemma [4. 3[ which gives 

P(AC T l > 0, AH Tu = 0) = 0, u > 0. (5.3) 

Let 

V H (dv) = / V(dt,dv)= / P(H s edv)ds 
Jt>o Jo 

be the potential measure of if, and similarly for Vg. 

Proof of Theorem 13.41 Assume X is not compound Poisson. Then is diffuse on 
(0, oo) by Lemma 1 of Chaumont and Doney (2010). Then since Hx({0}) = we have 

/ / V(ds,dv)U x ({v})= [ V s (dv)U x ({v}) = 0. (5.4) 

Js>0 Jv>0 Jv>0 
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Now by dSHUD , 

U. L -\H(ds, dx) = / V(ds, dv)U x (dx + v) for all s > 0, x > (5.5) 

Jv>0 

is equivalent to 

n. L -i tH {ds, {0}) = / V(ds, dv)U x {{v}) for all s > (5.6) 

Jv>0 

which in turn, by (15.4p . is equivalent to n L -i ,#(ds, {0}) being the zero measure. This is 
equivalent to X not creeping by Theorem 13.31 

If X is compound Poisson, then recalling the definition of (L~ l ,H) prior to (12. 6p in 
this case, a natural approach is to apply the result in the non compound Poisson case to 
the approximating process X\ = X t + et and take the limit as e I 0. Unfortunately this 
cannot work since X £ creeps and so the non compound Poisson result does not apply. 
Consequently we are forced to appeal to a direct construction of (L -1 ,if) in this case. 
We defer the details of this to the appendix. □ 

Example 5.3 Assume that X is a spectrally negative compound Poisson process with 
positive drift dx- For simplicity, to avoid killing, also assume EX\ > 0. In this example 
X creeps, and it is a simple matter to find IT^-i ,jj-(ds, {0}). Choose the normalisation of 
local time so that 

L t = J l(X r = X T )dr. 
Jo 

Let fTi be the time of the first jump of X and 

a = inf{t ><7i :X t >X ai }. (5.7) 
Then it is easy to see that (draw a picture) 

L~ x = s + y Rj, H s = d x s, 

i=l 

where iV is a Poisson process of rate A = Ilx(K) independent of the i.i.d. sequence 
Ri,i > 1, where R\ = a — o\. From this we conclude 

n L -i iH (ds, {0}) = \P{R X E ds) = U L -i(ds), s > 0. 

On the other hand, since X is not compound Poisson, (15. 4p shows that the lefthand side 
of (I3.10P is the zero measure when x = 0. 



Proof of Theorem 13. 5t Upon using ( I4.2p . this is an immediate application of Theorem 

S31 □ 
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6 Appendix; Proof of Theorem 13.41 in the Compound 
Poisson Case 



We want to prove Theorem 13.41 in the compound Poisson case. As mentioned earlier the 
natural approach using the approximating process Xf = X t + et does not work. Thus we 
proceed by using the random walk embedded in X to give a direct construction of the 
bivariate ladder processes. 

Throughout this section we assume that X is compound Poisson. We write it in 
the form X t = Ylk=i ^k, where N is a Poisson process of rate A > 0, Yh are i.i.d. rvs 
independent of iV with distribution function F, and P(Yk = 0) = 0. Thus the Levy 
measure of X is ITs^dx) = AF(dx). 

Let So = 0, and S n = Y^l=i ^k, n > 1. Let cr = and cr^, k > 1, denote the successive 
jump times of X t . Then {S n ,n > 0} is independent of {<Jk,k > 0} and X t = S n for 
o~ n < t < cr n+1 , n > 0. Let t = 0, 

t n+1 = mm{m > t n : S m > S tn }, n = 0,1,2,..., (6.1) 

be the weak ascending ladder times of S n , h n = S t „ the corresponding ladder height 
sequence, and 

U(k, x) = p ( f n = k,h n <x), k > 0, x > 0, 

n>0 

the corresponding bivariate renewal measure. The strict ascending ladder process is ob- 
tained by replacing the inequality in ( 16 .ip with a strict inequality. The analogous quan- 
tities for the strict descending ladder process will be denoted with a hat; thus 

U(k,x) = ^P(f n = k,hn <x), k > 0,x > 0. 

n>0 

We choose the normalisation of the local time L so that 

L t = / l(X r = X r )dr= / ^l Kii(Ttn+l) (r)dr, t>0. (6.2) 
</o Jo n=0 

As remarked earlier this gives rise to the weak bivariate ladder process (£ _1 , H) of X. For 
X, we require (L -1 ,if) to be the strict bivariate ladder process, which may be viewed 
as the limit of the ascending bivariate ladder process of X t — et as e 4 0. For a direct 
construction of (L~ l , H), let M s be an auxiliary Poisson process of rate 1, independent of 
X, and set 

jj \ _ i{ a t n ^n) if M 8 = n,t n < 00, n > 0, 
1 (00, 00) if M s = n, t n = 00, n > 0. 

This is analogous to the definition of the ascending ladder processes of X in the non- 
compound Poisson case when is irregular for [0, 00) for X; see for example p. 24 of 
Doney (2005). To emphasize, V(-, •) and •) are defined in terms of (£ _1 , H) given by 
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(16.31) . We should also remark that with these normalisations of the local times, one can 
check that (J2l| holds. 

The connection between the bivariate ladder processes of X and £ is given by 

Lemma 6.1 For any t > 0, x > 0, 

V(t, x) = A" 1 p (^ + i < t)U(k, x), (6 4) 

fc>0 

and 

V(t,x) = J2 P (°k<t)U(k,x). (6.5) 

fc>0 

Proof of Lemma I6.lt Since (16 .4p will not be used in the sequel, we only prove (16.51) . 
We have 

oo 

P(L; 1 < t, H s < x) = P(*u <t,h n < x, M s = n) 

n=0 

oo oo 

= ^2^P(t n = kJi n <x,a k <t,M s = n) (6.6) 

n=0 fc=0 

oo oo 

= P (*» = kX< x)P{cx k < t)P{M s 



— — n ) 

k=0 n=0 



One easily checks that J °° P(M S — n) ds = 1 for every n > 0, hence integrating (16. 6p we 
obtain 



V(t, X ) = J2J2 P & = k ^n< X ) P fa < t) 
k=0 n=0 

= J2P((7 k <t)U(k,x), 



k>0 

as required. □ 
Recall the definition of a in (15.71) . which in the present case reduces to 

a = inf{t >ai : X t > 0}. 

Also introduce 

(3 = M{n > : S n > 0}. 
Lemma 6.2 For s > 0, x > and v > 

V(ds, dv)F(dx + v) = P(X Q _ G -dv, X a G dx, a - a x G ds). (6.7) 
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Proof of Lemma 16. 21 Using P(<Jq G ds) = when s > in the fourth equality below, 
duality in the fifth, and (16. 5p in the sixth, we have 

P(X a ^ G — dv ,X a G da;, a — <J\ G ds) 

= l{ v>0 }P(X a _ G — dv, X a G dx, a — o x G ds) 

= 1{. U>0 } 2j P{P = h Si-i G -dv, Si G dx, Oi- cti G ds) 
= 1 { , >0} ^ P(/9 > i - 1, SU G -dt;)F(dx + v)P(a l _ 1 G ds) 
= l {u >o} Yl P (P > ^ e -dv)P(<n G ds)F(dx + v) 
= l {v>0} & (h dv)P(<Ti G ds)F(dx + v) 

i>0 

= l{. u>0 }U(ds, dv)F(dx + v) 
= V(ds,dv)F(dx + v) 

since s > and V(ds, {0}) = 8 (ds) by (E5]) (or (Q). □ 

Proof of Theorem 13. 4L Compound Poisson case: Since X does not creep in this 
case, we need to prove (13 . lOj) for all s > 0,x > 0. When s = x = 0, (13.101) reduces to 
showing 

/ V({0},dv)U x ({v})=0. (6.8) 

Jv>0 

But y({0}, df ) = 5o(dv ) by (16. 5p . which proves (16. 8p . Thus we may assume s V x > 0. 
If s > and x > 0, then integrating out t> in (16. 7\i gives 



y(ds, dv)F(dx + v) = P(AL; i i G ds, AH ai G dx). (6.9) 

) 

This continues to hold when s = and x > since 



V({0}, dv)F(dx + v)= 5 (dv)F(dx + v) = F(dx), 

v>0 Jv>0 

while 

P(AL~l = 0, AH ai G dx) = P(X ai G dx) = F(dx). 
Thus (16. 9 P holds whenever s V x > 0. Now by the compensation formula (which requires 
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s V x > 0) we have 

P(AL; i 1 G ds, AH ai edx) = Ej2 H L I- = t, H t . = 0, AL; 1 G ds, AH t G dx) 

= / P(L^ 1 = t,H t = 0)dt U L -i tH {ds, dx) 
Jt>o 

= / P(a 1 > t) dt U L -i >H (ds, dx) 
Jt>o 

= A _1 n L -i ijH '(ds,dx), 

which together with Ux(dx) = XF(dx) completes the proof of ( 13.101) . □ 
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